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Abstract. In this paper we construct a parametrix for the forward fundamental solution 
of the wave and Klein- Gordon equations on asymptotically de Sitter spaces without caustics. 
We use this parametrix to obtain asymptotic expansions for solutions of (□ — X)u ~ f and 
to obtain a uniform estimate for a family of bump functions traveling to infinity. 



1. Introduction 



De Sitter space is an exact solution of the vacuum Einstein equations with positive cosmo- 
logical constant. In this paper, we study the forward fundamental solution of the wave and 
Klein-Gordon equations on asymptotically de Sitter spaces. This is the unique operator 
(which we identify with its Schwartz kernel) which satisfies (□ — X)E^ = I and is supported 
in the forward light cones, i.e., for a compactly supported smooth function /, the function 
u = E^f satisfies 



Here A is the Klein-Gordon parameter. If equation ([T]) is considered as a massive wave 
equation, the condition A < corresponds to positive mass. We construct a parametrix 
for this problem and establish asymptotic expansions for solutions. As an application of 
our parametrix, we prove a uniform estimate for the operator applied to a family of 
bump functions tending toward infinity. We postpone to a future paper the consideration of 
Strichartz estimates and semilinear wave equations on asymptotically de Sitter spaces. 

The study of the decay properties of the wave equation on various natural classes of 
spacetimes is an active area of research. For example, Dafermos and Rodnianski [5] and 
Melrose, Sa Barreto, and Vasy [17] have obtained decay results for solutions of the wave 
equation in the context of the de Sitter-Schwarzschild model of a black hole spacetime. 

Our definition of asymptotically de Sitter spaces is given in [21] (called asymptotically de 
Sitter-like spaces in that manuscript) and follows the definition of asymptotically hyperbolic 
spaces given in [T^ and [16]. An asymptotically de Sitter space is a compact manifold with 
boundary equipped with a Lorentzian metric having a prescribed asymptotic form near the 
boundary. This pushes the boundary off "to infinity" . 

The microlocal structure of the fundamental solution for general real principal type oper- 
ators has been studied extensively. The solution operator for the Cauchy problem for gen- 
eral real principal type operators is a Fourier integral operator associated to a Lagrangian 
submanifold of phase space given by the fiowout of the Hamilton vector field of the prin- 
cipal symbol of the operator. This was first described in this language by Duistermaat 
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and Hormander in [5]. In [18], Melrose and Uhlmann constructed the forward fundamental 
solution for a real principal type operator as a paired Lagrangian distribution. These are 
distributions associated to two cleanly intersecting Lagrangian submanifolds in phase space. 

Guillemin and Uhlmann [5], Joshi [TT], Melrose and Zworski [TH], Hassell and Vasy [9], 
and others have all generalized the notion of paired Lagrangian distributions. Guillemin and 
Uhlmann defined a much more general class of paired Lagrangian distributions, which Joshi 
restricted slightly in order to construct a well-behaved calculus. Joshi then used this calculus 
to construct complex powers of the wave operator on Riemannian manifolds. Melrose and 
Zworski defined a class of distributions associated to an intersecting pair of Legendrians, 
while Hassell and Vasy later expanded this notion to describe the spectral projections on a 
scattering manifold. 

Polarski |_23j computed the propagator for the equation of the massless conformally coupled 
scalar field in the static de Sitter metric, which has been transformed to the Einstein open 
universe. Yagdjian and Galstian [25j computed the fundamental solutions for the Klein- 
Gordon equation in de Sitter spacetime transformed by the Lemaitre-Robertson change 
of coordinates to the special case of the Friedmann-Robertson-Walker-Lemaitre spacetime. 
They represented the fundamental solutions and solutions of the Cauchy problem using 
hypergeometric functions and proved VL'^ estimates. 

Vasy [21] generalized and extended Polarski's result to asymptotically de Sitter spaces. 
He exhibited the well-posedness of the Cauchy problem and showed that on such spaces, 
the solution m of (□ — X)u = with smooth Cauchy data has an asymptotic expansion at 
infinity. Indeed, if x is a boundary defining function for the conformal compactification of 

an asymptotically de Sitter space X and \J~^^-^rX is not a half-integer, then u has an 
expansion 

U = U+X 2 ^ V 4 + U-X 2 V 4 ^^^^ 

where u+ and U- are smooth on X. (The difference between this expression and the corre- 
sponding one in [21] is due to differing sign conventions for Lorentzian metrics.) In the case 

/ („_i)2 

of an integer coincidence, U- is instead in C^{X) + x^V ^^+^(loga;)C~(X). Vasy also 
showed that solutions of the wave equation exhibit scattering, i.e., that the data u± may be 
specified at one of Y± (future and past infinity, respectively), which fixes the data at Y^. 

Our result extends the work of Vasy to the study of the fundamental solution (□— A)i?+ = / 
on asymptotically de Sitter spaces. We require three global assumptions in our study of 
asymptotically de Sitter spaces: 

(Al) y = y+ U y_, with Yj^ and Y_ a union of connected components of Y , 
(A2) each bicharacteristic 7 of P converges to y+ as t — )■ +00 and to y_ as t — )■ —00, or 
vice versa, and 

(A3) The projection from T*X to X of the flowout of the forward light cone from any point 
p G YL is an embedded submanifold of X (except at the point p, where it always has 
a conic singularity). 

The first two of these assumptions are not particularly restrictive. They imply that the 
manifold is topologically a product Y^ x M, but are reasonable from a physical viewpoint 
in that they imply a time orientation on the manifold. Colloquially, assumptions (Al) and 
(A2) prevent the breakdown of causality on X. 
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Assumption (A3) is needed only to obtain sharp global statements about the fundamental 
solution, but our construction works in some neighborhood of (i.e., a neighborhood of 
future infinity) without this assumption. 

Assumption (A3) prohibits the development of caustics, which significantly narrows the 
class of manifolds considered. De Sitter space just misses being covered by (A3), though a 
slight modification of our construction still applies here. Indeed, the projection of the flowout 
of the light cone from a point p G F_ is a smooth embedded submanifold of the interior of de 
Sitter space, but intersects itself at Y+. Section [T6] discusses the minor modifications needed 
to handle this case. 

If we slightly enlarge the spherical cross section of de Sitter space, then this new space 
satisfies the assumptions above. The assumptions above are stable under perturbation, so 
the construction applies to perturbations of this enlarged version of de Sitter space. 

In order to remove the assumption (A3), we could combine the Poisson operator construc- 
tion of [21] with the local description of the forward fundamental solution (given in |18j). 
Phrasing this in a geometric way is left to a future paper. 

The main result of this paper is the following (we state it more precisely and define the 
relevant classes of distributions later): 

Theorem 1. Suppose {X,g) is an asymptotically de Sitter space. There is a compactification 
Xq of the interior of X x X to a compact manifold with corners such that the closures 
of the diagonal and the light cone both intersect all boundary hypersurfaces transversely. 
Xq is constructed by first blowing up the boundary of the diagonal in X x X and then 
blowing up the set where the projection of the flowout of the light cone hits the side face. 
On this compactification, the forward fundamental solution o/ □ — A lifts to be the sum of a 
paired Lagrangian distribution, smooth down to the front face, and a conormal distribution 
associated to the light cone with polyhomogeneous expansions at the other faces of Xq. 

We call a tempered distribution f on X forward- directed if it a smooth function on the 
interior of X that vanishes to all orders at Y_. The work of Vasy [23] implies that if / is 
forward-directed, then we may apply to /. 

A function / is polyhomogeneous with index set £" on X if it has an asymptotic expansion 
of the form 



A corollary of Theorem [T] describes the polyhomogeneity of the solutions of P{X)u = f. 

Theorem 2. // / is forward- directed and polyhomogeneous on X with index set E at Y+, 
then Ej^f is forward- directed and polyhomogeneous with index set F , where 




(r,l)GE 





F = F1UF2UE. 
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Here is regarded as an operator and U denotes the extended union of two index sets: 

{z,p) G EOF ^ {z,p) e E U F OT p = p' + p" + 1 with {z,p') G E and {z,p") G F. 

As an apphcation of our parametrix, we obtain a uniform estimate for a family of bump 
functions travehng to infinity. We leave more general mapping properties and Strichartz 
estimates to a future manuscript. 

Theorem 3. Suppose (p G C^(]R^ x M""^) is supported near (1,0). For {x,y) G X, Let 
f{x,y){x,y) = Suppose that p, I, and r satisfy 

2 < p < oo, 




r > max | h A | , and 



1 n-l ' 
2' 2 , 



I > max 0, — + - — ^-^ + A 

Then 

(3) \\^+fi^'y)\\x-l+^('~-)/PLP{X;dg) - 

Here the constant C depends on the family f(x,y) hut not on {x,y). 

In particular, for the wave equation, A = 0, and we may let I = and r > max ( 
In this case, 

II F f,- -^ II < r^~'^'^lp 

Theorem [1] follows from a parametrix construction that comprises the bulk of this paper. 
The main ingredient in the proof of Theorem [2] is the pushforward theorem of Melrose. 
Theorem in] follows from Theorem [T] and the estimates proved by Vasy |24j . 

The paper is broadly divided into three parts. The first part consists of Sections [2] through 
[3 and develops the tools necessary to construct the parametrix. The second part contains 
the construction, which begins with an outline in Section [TU] and concludes in Section [T3 
The last part, consisting of Sections [T71 and [T8| proves Theorems [2] and [31 



2. Asymptotically de Sitter spaces 

We start by describing the de Sitter space. Recall that hyperbolic space can be realized as 
one sheet of the two-sheeted hyperboloid in Minkowski space. It inherits a Riemannian metric 
from the Lorentzian metric in Minkowski space. De Sitter space, on the other hand, is the 
one-sheeted hyperboloid {— Xq + Y^^=i-^i = 1} in Minkowski space, but now the induced 
metric is Lorentzian. A good set of coordinates on this space, which is diffeomorphic to 
5" X M, is 

Xq = sinh r 
Xi = uji cosh r, 

where Ui are coordinates on the unit sphere. The de Sitter metric is then 

— dr^ + cosh^ r duj"^. 
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If we let T = e ^ near r = +00, then r = +00 corresponds to T = and the metric now 
has the form 

-dT^ + \{T^ + lfduj^ 



J^2 

Our definition of an asymptotically de Sitter space is based on the form (j4]) of the de Sitter 
metric. 

Definition 4. {X, g) is asymptotically de Sitter if X is an n- dimensional compact manifold 
with boundary Y, g is a Lorentzian metric on the interior of X, and, for a boundary defining 
function x, there is some collar neighborhood [0, e)x x y of the boundary on which g has the 
form 

dx"^ h{x,y,dy) 
9 = 5- 2 ' 

where h is a smooth symmetric (0, 2)-tensor on Y that is a metric on the boundary {0} x Y. 



A calculation similar to the ones in [2] or [TB] shows that the sectional curvatures of 
asymptotically de Sitter spaces approach 1 as x — )■ 0. 

Let □, the wave operator, be the wave operator (d'Alembertian) associated to g, and, for 
A G M, let P = -P(A) = □ — A be the Klein-Gordon operator. We are seeking the fundamental 
solution of the Klein-Gordon equation, i.e., a distribution E+^\ such that P{X)E^^x = I. Note 
that the principal symbol (72 (P) is given by the dual metric function of g. 

Definition 5. We say that the bicharacteristics of P over X°, the interior of X, are the inte- 
gral curves of the Hamilton vector field of the principal symbol o'2{P) inside the characteristic 
set of P (the set where cr2(P) vanishes). 

Throughout this paper, we assume that both (Al) and (A2) hold. These assumptions are 
introduced primarily to ensure that X exhibits a coherent causal structure. 
Because g is conformal to the incomplete pseudo-Riemannian metric 

g = —dx^ + h 

near Y , the bicharacteristics of P (near Y) are reparametrizations of the bicharacteristics of 
—dx'^ + h. g is complete, so the global assumptions imply that each bicharacteristic 7(t) has 
a limit in SyX as t — )■ ±00. 

Physicists (e.g. [B]) have long known that these assumptions imply the existence of a global 
'time' function T G C°°{X) such that T|y^ = ±1 and is monotone on the nuUbicharacteristics 
of P. Note that 1 — x and x — 1 have the desired properties near y+ and Y^, so the 
assumptions mean that such functions can be extended to all of X. Moreover, T gives a 
fibration X — )■ [—1, 1] and so X is diffeomorphic to X x for a compact manifold 5*. This 
also shows that Y^ and Y_ are both diffeomorphic to S. 

We also assume that assumption (A3) holds. This assumption ensures that the boundary 
of the flowout of the light cone in the cotangent bundle is actually the cotangent bundle 
of a submanifold in the base. It allows us to work with an adapted class of conormal 
distributions rather than a class of Lagrangian distributions. Note that the projection of the 
flowout is automatically an embedded submanifold for small times. Moreover, because the 
nuUbicharacteristics agree in the compact and non-compact settings (i.e., for g and g), there 
is a neighborhood of Y^ such that the projection of the flowout of the forward light cone 
from any point in this neighborhood is an embedded submanifold. If we then restrict to data 
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supported in this neighborhood, we may use our construction below even in the absence of 
assumption (A3). 

In this paper we adopt the convention that / 5{x — x)6{y — y)f{x,y)dg = f{x,y). In 
particular, this breaks the formal self-adjointness of the operator -P(A). When we seek a 
right parametrix, then, we must use the transpose of the operator P(A) with respect to this 
metric. We may write this out explicitly as 

(5) P\X) = (xd,)' + (n + l)xd, + ^-^^xd, + x'^y + n{\ + - A. 

We also adopt the convention that when we are applying a differential operator to the 
variables in the right factor, we use a subscript R. 

In our consideration of the forward fundamental solution of -P(A), it is useful to have the 
notion of forward- directed data, which are simply tempered distributions (or polynomially 
bounded smooth functions) on X that vanish to all orders at Y_. 



3. The Cauchy problem and scattering 

In [21], Vasy studied both the Cauchy problem and the scattering problem on asymptoti- 
cally de Sitter spaces. The Cauchy problem seeks a solution the equation 



(6) 



|p(A)w = 0, 



where Sq is a Cauchy hypersurface, V is a vector field transverse to Sq, and ipo,4'i are 
smooth functions on Sq. 

Via a positive commutator argument, Vasy showed global existence and uniqueness for 
solutions of the Cauchy problem ([6]) on asymptotically de Sitter spaces. One may use any 
"time slice" {T = const} from the diffeomorphism X = M x "K^^ as the Cauchy hypersurface 
needed to pose the Cauchy problem. 

Vasy further proved the following a global solvability result for the inhomogeneous equation 
on asymptotically de Sitter spaces under assumptions (Al) and (A2). 

Theorem 6 (Theorem 5.4 of |2l])- Suppose that A G M, and that l± satisfy 

/+ > max(-, /(A)), /_ < — max(-, /(A)). 
Then for f E Hq'^'^'''' (X) , Pu = f has a unique solution u E Hq^^'''^ (X) . 



Here /(A) = 'Siy + A while the if^'''+''' = x'^_^^x'L measure both regularity and 

decay at separately. The spaces are weighted Sobolev spaces measuring regularity 

with respect to the 0-vector fields defined in the following section. 

In particular, for any / E C°°{X) (indeed, for any forward-directed tempered smooth 
function), there is a unique solution u E x^^x'^Hq{X) to Pu = f. (Here x± is a defining 
function for Y±, so the notation means that u is tempered at and vanishes to infinite 
order at Yl.) There is thus a distribution on X x X that can be called the forward 
fundamental solution of P. In other words, is such that for each p E X, PE^{p) = 6p 
and E+{p){q) = for q E X not in the domain of influence of p. In terms of the identification 
X = X F+, this means that E+{p){q) vanishes when t{q) < t{p). 
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Vasy (cf. Theorem 5.5 of [24j) proved also that solutions of the Cauchy problem ([6]) have 
unique asymptotic expansions near the boundary F of X of the form 

U = U + X 2 V 4 -t--^ + 2 V 4 "t--^^ 

where m|± are smooth up to the boundary of X. When the difference between the two expo- 

nents is an integer, then U- is instead only in C°°{X) + x V 4 + (\ogx)C°°{X). Moreover, 
one may specify u±\y_, uniquely determining u±\y^. He showed (cf. Theorem 7.21 of 
that when there are no integer coincidences the map sending the data 

{u+\y_,U-\y_) {u+\y+,U-\y+) 

is a Fourier integral operator associated to the canonical relation given by bicharacteristic 
flow from YL to Y^. This is via an explicit parametrix construction on an appropriate blow- 
up of X x Y, and the static model of de Sitter space appears in the consideration of the 
rescaled normal operator at the front face of this blow-up. 

Observe now that for [xq, yo) G X away from the boundary of X, the restriction of the 
kernel of the fundamental solution to the slice x = Xq agrees with a multiple of the solution 
of the Cauchy problem 

P{\)u = 0, 

u{xo, ■) = 0, 

Vu{xo, ■) = 5{y - y), 

at least in the region to the future of (xo,yo)- (Here V is as in equation ([6]), i.e., a linear 
combination of xdx and xdy transverse io {x = xq}.) In particular, we can understand the 
fundamental solution in this region by understanding the solution operator for the Cauchy 
problem, which Vasy studied. 

The solution operator for the Cauchy problem is the composition of the Poisson operator 
and the Cauchy-to-scattering operator. Vasy showed that these are Fourier integral oper- 
ators, and their canonical relations intersect transversely, so we can compose them. The 
result is an operator with canonical relation given by the restriction of the conormal bundle 
of the flowout of the light cone, restricted to a; = a;o but lifted to a blown-up space. This 
blown-up space agrees with a slice of the space we define later. Vasy's construction blows 
up [X X y+,diagy+], which turns into what we call lcf+ under this composition. Applying 
this operator to a delta distribution gives a conormal distribution on our space, and the log 
terms in Vasy's construction become the log terms in our construction, though this requires 
careful bookkeeping. Viewing this as a distribution parametrized by x gives the fundamental 
solution in a neighborhood of the interior of what we call lf_(. away from the other boundary 
hypersurfaces. 

4. 0-GEOMETRY 

Recall from [16] the Lie algebra of 0-vector fields, 

Vo = {vector fields vanishing at dX}. 
The universal enveloping algebra of this Lie algebra is the algebra of 0-differential operators, 

Difr*(x). 
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We may use the Lie algebra Vo{X) to define the natural tensor bundle '^TX, whose smooth 
sections are elements of Vo{X). Its dual, ^T*X, has smooth sections spanned (over C°°{X)) 
by ^ and ^. 

Recall that T*X is endowed with a canonical 1-form given by 

/- , , ^dx dy 
a = t,dx + jjL ■ dy = xt, h x/i ■ — . 

X X 

In particular, ^T*X is endowed with the canonical 1-form 
(7) o^ = r— + r]-^. 

X X 

Just as T*X is endowed with a symplectic form given by u; = da, '^T*X is endowed with 
the symplectic form = d^a. 

As in [T3] or [16], we define the 0-double space X Xq X = Xq as the blown- up space 
[X X X, dA]. X Xq X is a manifold with corners agreeing with X x X on the interior. It has 
three boundary hypersurfaces - the left face If = {x = 0}, which is the lift of the boundary 
hypersurface {0} x X in X x X; the right face rf = {x = 0}, which is the lift of X x {0}; and 
the front face ff, which is the boundary hypersurface introduced by the blow up construction. 
Recall also that Xq is equipped with a blow-down map /3 : Xq — )■ X x X. A neighborhood 
of the front face ff is depicted in Figure [H 




Figure 1. Passing from X x X to the 0-double space Xq = [X^, dA]. 

This construction is perhaps best thought of as an invariant way of introducing spherical 
coordinates near the boundary of the diagonal in X x X. Indeed, a valid coordinate system 
near the front face ff is given by spherical coordinates: 



Pff 



a/x^ + x^ + ly — 



n _ ^ . _ 5 n> _ y-y' 

tfo — , t/n — , O — , 

Pff Pff Pff 

where 6 G 5*"+, a quarter sphere. It is often more convenient to work with projective 
coordinates near the front face: 

_ X _ y-y ^ , 
s _ , z _ , x,y. 

X X 

In the projective coordinates {s,z,x,y), we may compute the lifts of the vector fields 



xdx, xdy, xdx, and xdy. Indeed, we have 



xdx ~^ sds xdy -w sdz 



xdx xdx — sds — z ■ dz xdy xdy — d^- 
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Recall that the front face is the total space of a fibration. Indeed, it is the total space of a 
bundle over Y with quarter-sphere fibers. In local projective coordinates (s, z, x, y) near the 
front face, the fibers are given by {x = 0,y = const}. An asymptotically de Sitter metric 
on X induces a Lorentzian metric — + on the interior of the fibers. By an affine 

change of coordinates in the fiber, this metric may be written as — + which is the 
Wick rotation of the (negative definite) hyperbolic metric on an upper half-space. 

Due to our assumptions (Al) and (A2), the space X has two boundary components F+ 
and y_, so it is useful at this point to give names to the components of the left face, right 
face, and front face. We call ff+ the component of the front face that comes from blowing 
up the diagonal near x Y^, and ff_ the other component. Similarly, we give the name lf+ 
to the lift of the hypersurface y+ x X and lf_ to the lift of y_ x X (and rf_|_, rf_ denote the 
lifts of X X y+ and X x YL, respectively). This naming scheme is illustrated in Figures [2] 
and 131 




Figure 2. A two-dimensional view of Xo. 




Figure 3. Another view of Xq. The horizontal lines represent the y,y axes. 

We should also recall the trivial half-density bundles ^^(X) and ^^(Xq), which are trivial 
bundles. In local coordinates {x,y), Q2(^X) is trivialized by the global section ■y = \dg\ ' . 

J 1 /2 

The bundle f22(X x X) is trivialized by t; = | dgidgjil , where dg^ and dgji are the lifts 
of the densities from the left and right factors of X x X. Up to a nonvanishing factor, the 
Jacobian determinant of the blow-down map Xq — )■ X x X is r", so t; lifts to r"/^/i, where 
is a nonvanishing section of the standard half-density bundle ^^(Xq). For example, we 
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~ 1 /2 

may take /i = \ drd9dy\ in local polar coordinates. In projective coordinates, we may 
take = \ds dz dx dy\^^'^ , and then v lifts to (times a nonvanishing factor). In this 

paper, we adopt the convention that v and 7 are flat, i.e., for a differential operator L, 
L{u'^) = {Lu)'~f and L{uv) = {Lu)v. In particular. 

Our operator P(A) commutes with x, so P{X){ufi) = {P{\)u)fi. 

On the fibers of the front face, elements of DiffQ(X) lift to differential operators on Xq 
that restrict to differential operators on the fibers of the front face by setting x = 0. We 
call this restriction to the fiber over p ^ Y the normal operator and denote it as Np{A) for 
A G Diffg(X). As an example, in projective coordinates the normal operator of -P(A) is given 
by 

(8) iVp(P(A)) = (sdsY - (n - l)isds) + s^A, - A. 

The inclusions ff C Xq and If C Xq induce inclusions on the tangent bundles T ff C TgX^ 
and Tlf C T\{Xq. These inclusions induce natural projections on the cotangent bundle: 

vTff : T;X^ ^T*S and 7TiiT*X^ T* If . 

Let 2Xg denote the doubling of Xq across ff , and let TqXq denote the restriction of the 
cotangent space to the front face. As in [12], we say that a smooth closed conic Lagrangian 
submanifold A C T*Xq is extendible if it intersects T^Xq transversely. Recall that in this 
case there is a smooth closed conic Lagrangian submanifold A*^ C T*2Xq such that 

A = A"nT*X2, A° = 7rff(AnTf;Xo^). 

Note that there are many choices for the extension A^. 
We recall a result [12] : 



Lemma 7 (Lemma 2.1 of [12]). Let A C T*Xq he an extendible Lagrangian. Then Aq = 
7rfj(A n T^Xq) is a Lagrangian submanifold ofT* ff . 

For future use, we collect a bit more information about the symplectic structure of T*Xq. 
The space ^T*X x ^T*X is endowed with the symplectic form 

(9) u = Tclux + rc*ujx, 

where Ux is the symplectic form on ^T*X (coming from the canonical 1-form ([7])) and ttj 
the projection of the jth copy of °T*X. Moreover, the blow-down map /3 induces a smooth 
map 

°T*X X °T*X ^ T*X^ 

which is an isomorphism on the interior. 

The identification of T*X and '^T*X over the interior of X then induces a smooth map 
^T*X X or*X — )■ T*Xq over the interior of X x X. This map identifies the 1-forms 

dx _ dx dy _ dy , ^ ,~ ^ , ,~ 

T h r— ^ + rj ■ h rj ■ and a ds + t, dx + Q ■ dz + fi ■ dy, 

•Xj ju ju 

and so we must have 



(10) 


a = 


TX 
X 


e = 


X X X X 




c = 


XT] 
X 


/i = 


1 + 1 

X X 
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Observe that in these coordinates, the symplectic form oj takes a famihar form: 

da /\ ds + dC, A dx + d( A dz + dfi A dy. 

5. POLYHOMOGENEITY AND CONORMAL DISTRIBUTIONS 

In order to consider the pushforward of a distribution, we need the notion of a b-fibration. 
We recall from ^20j that an interior h-map f : X oi manifolds with corners is a function 
mapping X — )■ F and X\dX — )■ Y\dY such that each boundary defining function for Y pulls 
back to a sum of products of boundary defining functions for X. More precisely, suppose 
that Mi{X) is the set of boundary hypersurfaces of X, and, for each G G Mi(X), let I{G) 
be the ideal of functions vanishing at G. Suppose that if is a boundary hypersurface for Y 
and I{H) is the ideal of functions vanishing on H. We say that / : X — )• y is an interior 
b-map if 

f*i{H) = W i{Gyf^^'"\ 

Here ef{G, H) is a collection of nonnegative integers that we call the exponent matrix of /. 
We set (e/) = {H e Mi{X) : e/(/7, G) = OVG E Mi{Y)}. 

For an interior b-map /, the differential /* : TpX — )■ Tf{p)Y extends by continuity to the 
b-differential 

: ^TpX -)■ ^Tf^p-)Y. 

Recall that a b-fibration between two manifolds with corners is one that is both b-normal 
and a b- submersion. A b-submersion is a b-map with surjective b-differential. A b-normal 
map is one such that the b-differential is surjective as a map ^N^H — )■ ^Xj(^)G on the interior 
of each boundary hypersurface H . Here G is the face that contains the image of H . 

We say that a discrete subset _E C C x No is an index set if {zj, kj) G E with \ {zj, kj) \ — )■ oo 
implies that ^zj — )■ oo, [z, k) E E implies that (2; -|- p, fc) G E for all p G No, and if 
{z,k) G E, then G E for all p G Nq with {) < p < k. We say that £ = {Eh : 

if is a boundary hypersurface of M} is an index family if each E^ is an index set. 

Our construction below shows that certain distributions have asymptotic expansions. We 
make this notion more precise by recalling the definition of a polyhomogeneous distribution. 

Definition 8. A polyhomogeneous distribution is a distribution that is smooth on the interior 
of M, and, near each boundary hypersurface H oi M has an asymptotic expansion of the 
form 

^ x^C^og XHyasi{x,y), 

where is a defining function for H, Eh is an index set for H, and a^z are smooth functions 
independent of xh- Near the corners of M, we require that polyhomogeneous distributions 
have an appropriate product-type expansion. 

For an extended discussion of polyhomogeneous distributions, we refer the reader to [15], 
[20] , or [21] . We adopt the convention that a^; should take values in the half-density bundles. 

Recall that if / and g are polyhomogeneous with index families J-" and Q, then fg is 
polyhomogeneous with index set J-' + Q. 

We also require the following lemma, which can be found in [1U\ or [2T], and allows us to 
understand the polyhomogeneity of the pushforward of a distribution: 
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Lemma 9. If f : X ^ Y is a b-fibration between compact manifolds with corners, then for 
any index family /C for X with 

^IC{H)>-1 ifHeicf), 

the pushforward gives 

/. : A^,^{X; n ® rE) ^ Q^E), 

J = f#ic. 

Here J7 = /#/C is the pushforward of the index family /C. In particular, for a boundary 
hyper surf ace H of Y, 

J{H) = OGeAMX),efiG,H)^o{{z/ef{G,H),p) : {z,p) G /C(G')}, 

where U is the extended union defined by {z,p) G KUJ if and only if {z,p) G K U J or 
p = p' + p" + I with {z,p') G K and {z,p") G J. 

We recall the definition of the space of conormal distributions (see, for example, or 
|2U]). If M is a manifold with corners and Z is an interior p-submanifold, then I^{M\ Z) is 
the space of conormal distributions of order p associated to Z. Here we only need the case 
when Z has codimension one and meets the boundary hypersurfaces transversely. In local 
coordinates (x, y) where M is given by Xj > for i = 1, . . . , k and Z is given by yi = 0, then 
u G P{M; Z) has the form 

e^^^^a{x, y, rj) drj, 

with a a half-density valued classical symbol of order p + n/A and n = dimM (note that 
we are using the order convention consistent with the orders of Lagrangian distributions in 
^18j). Because x behaves as a parameter, the boundaries cause no problems here. 

We also need several refinements of this notion. The first is when we allow the symbol 
a to have an asymptotic expansion (in x) near the boundary of M. We need the space 
ApYigI^{M] Z) of conormal distributions associated to Z with polyhomogeneous expansion 
at the boundary of M. These are merely conormal distributions in the sense above, but 
where we allow the symbol a to have a polyhomogeneous expansion in the x variables. 

The second generalization we require is to allow the symbol to have a mild type of singu- 
larity at 7] = 0. In particular, we allow symbols of the form Yljiv ~^'''^Y~''(^j- The singularity 
at only affects the growth of the distribution as ?/ — )■ ±oo, and so we may avoid compli- 
cations by multiplying our distribution by a cutoff function supported near y = 0. We need 
this generalization later in order to prove a lemma about the support of our distributions. 

6. The Hamilton vector fields 

By analogy with the definition of the Hamilton vector field, we recall from [T2] the defini- 
tion of the 0-Hamilton vector field. 

Definition 10. Given p G C°°(T*X), the 0-Hamilton vector field of p, denoted ^Hp, is 
defined by 

^ui;^Hp) = dp. 

In local coordinates where °a is given by ([7]), ^Hp is given by 

0^ dp d dp d f dp dp \ d f dp dp \ d 
^ dr dx df] dy \ dx dt] ^ / dr \ dy dr ^ J dr] 
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We are interested in the operator -P(A), whose principal symbol is the length function 

P = -x^i'^ + x'^ihoiy, jj) + xhi{x, y, /i)) = -r^ + ho{y, t]) + xhi{x, y, r]), 

where we have written h{x,y,fi) = ho{y,fi) + xhi{x,y, This is also the principal symbol 
of P when P is treated as an operator acting on the left factor of the product space X x X. 
In particular, in the coordinates given by (ITOl) . p pulls back to 

p = —s^a'^ + s'^ho{y + xz, () + s^xhi{xs, y + xz, (). 

Thus if u is the lift of the symplectic form (Q (as in [12]) to T*Xq, ^Hp lifts to Hp, which 
is given by a)(-, Hp) = dp. In these coordinates, we may write 

^ dp d dp d dp d dp d 

^ da ds ds da d^ dz dz d( 

dp d dp d dp d dp d 

d^ dx dx d^ dfi dy dy dfi 

An elementary computation then yields 



d I dh \ d 

Hji = - 2sV— - s -2^2 + 2ho + Ssxhi + s'^x^—^ — 
ds \ dx J da 

_^ ^2 f^j^ sx— ^ ■ — - (x— + sx^—^ ■ — 
\dr] df] J dz \ dy dy J d( 

2/ ^^0 , , , 2~9hi _ dhi\ d 2 f 9ho , ^pdhA d 
[z - —- + shi + s X— + sxz ■—-]—- s [^— + sx—— ■ — , 
\ dy dx dy J dt, \ dy dy J dfi 

where ho and hi are evaluated at {xs, y + xz, (). 

Note that the characteristic set of P on Xq is the set where p = 0, i.e., the set where 
s^cr^ = s^ihijj + xz, C) + sxhi{sx, y + xz, ()). In this set. Hp has the form 

2/ o f-i. ~2^^l^ "5 , f 9ho ^ ^dhi\ d 



Hp=s [ - 2a— - [xhi + sx — — ^ + ^— + sx , 

' ds \ dx J da \ dr] drj / dz 

_dhQ _2^^i\ ^ f ^^0 ^pdhi\ d 
— I X— h sx ^— ■ tt: — -T^ 1- sx- ' 



dy dy J dQ \ dy dy J dfx 

dho ^ , ^ 2~'9^i , ~ dhi\ d 



(11) - [ z ■ —- + shi + s X—- + sxz ■ 

\ dy dx dy J dt, 

In particular, we have proved the following proposition: 

Proposition 11. Inside the characteristic set of P on X^, Hi = Hp = s'^Hl, where Hl 
is a smooth vector field tangent to the front face that is nondegenerate at s = 0. Similarly, 
in coordinates {x,y,s = s~^,z = s~^z), the Hamilton vector field for the lift from the right 
factor can be written Hr = s^Hr, where Hr is nondegenerate at s = 0. 



The fact that we may factor a power of from the Hamilton vector field is useful in the 
next section. 
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7. The Lagrangians 

The two Lagrangians that interest us are the hft of the diagonal and the flowout from the 
characteristic set of P within this conormal bundle by the Hamilton vector field H^. 

The coordinates on phase space given by fITU]) imply that the conormal bundle of the 
diagonal, 

{x = x,y = y,T = -f,7] = -fj}, 

lifts to 

Ao = {s = l,z = 0,^ = 0,/i = 0}. 
The Lagrangian submanifold Aq intersects T^Xq transversely at 

A° = {s = l,z = 0,e = 0,/i = 0,x = 0} = T^^nffff, 

and so Aq is extendible. 

We now set Ai to be the forward flowout of Hp from Aq H S(P). Because we may write 
Hp = s'^Hp, with Hp a smooth nondegenerate vector field tangent to ff_|_, Ai is a smooth 
submanifold of T*Xq, intersecting Tg Xq transversely. Ai is Lagrangian by general theory 
(see, for example, [1]). Thus Ai is an extendible Lagrangian near ff+. 

To see that Ai is extendible near ff_, we require the following proposition: 

Proposition 12. The Lagrangian Aq is invariant under the flow of Hi — Hr. Moreover, Ai 
is also the flowout by Hr of Aq fl S(P). 

Proof. A straightforward calculation shows that Hl — Hr preserves Aq. Indeed, though this 
is a nonzero vector field, the coefficients of ds, dz, d^, and 5^ all vanish at Aq. Note further 
that S(Pi) n Aq = S(Pr) n Aq, and so H^ — Hr preserves this set. 

In order to see that Ai is also the flowout of iiT/j, we observe that H^ and H^ must 
commute because they are lifts (from X x X) of commuting vector fields. In particular, 
because Hr — Hi is tangent to Aq fl S(Pl), it is also tangent to Ai. In particular, H^ is also 
tangent to Ai. Hji is not tangent to Aq, so its fiowout must also be Ai. □ 



An argument similar to the one above then shows that Ai is extendible near ff_ (with Hl 
replaced by Hji). (For what we use below, the argument used above suffices, as we only need 
that Ai is extendible in a neighborhood of the diagonal.) 

The projection of H^ in S(P) to T* ff_|_ is given by 

^ \ ds d( dz dy djj, 

so that A^ is given by the flowout of this vector fleld from Aq in T* ff+. 

Note that the assumption (A3) implies that the projection of Ai to the base Xq is a 
smooth submanifold. We call this submanifold LC = vrAi. This is a smooth codimension 
one submanifold. Note further that LC is the boundary of an open subset of Xq, which we 
call LC*"*. 
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8. Paired Lagrangian distributions 

We are interested in a calculus of paired Lagrangian distributions adapted to the 0- 
geometry. These are distributions with a model form 







e'^'-^-'^^+'^-^a{t, s,Y,x, y, ^, jj.) dt djj, 



where a is a symbol. 

We briefly recall the definition of a paired Lagrangian distribution from [18]. Suppose M 
is a manifold and Lq,Li are conic Lagrangian submanifolds of T*M intersecting cleanly in 
codimension 1. We say that u G I'^{M; Lq, Li, Q2{M)) if it is a distributional ^-density 
modeled microlocally on 

Jo Jr" 

where a is a symbol of order m + ^ — j. In the general case, we require the phase function 
to parametrize the Lagrangian pair (Lo,Li) and a to be a symbol of order m + | + , 
where N is the number of variables required to parametrize Li. In this model, the phase 
function parametrizes Lq = A^*{0}, while Li = {xi >0,x' = 0,C,i = 0}. 

We now recall what it means to parametrize a cleanly intersecting Lagrangian pair. Sup- 
pose that Lq,Li C T*M are two closed conic Lagrangian submanifolds intersecting cleanly 
in codimension 1 and that dLi = Li n Lq. Suppose : X \ {0} — )■ M is homogeneous 
of degree 1 in ^, and let 

C = {(x,0 : xGM,eeM^\{O},rf50(x,O = O}. 

We say (h is nondegenerate if ci^, ? = 1, . . . , A^, are linearly independent at any point in 
C . We say that parametrizes a single conic Lagrangian submanifold L if L is locally the 
image of the map C T*M given by {x,C,) {x,dx4>{x,C,))- A simple example is the case 
when L = N*{0}, in which case we may take (j){x,^) = x ■ ^. 

In the case of an intersecting Lagrangian pair, we say that is a nondegenerate parametriza- 
tion of (Lq, Li) near q G dLi if 

s, = 0o(a;, + s, 0, ^ e M", 

where q = {xo,dx4>{x,0,^)), d(^s,^)(p{x,0,^) = 0, 0o is a nondegenerate parametrization of Lq 
near q, and parametrizes Li near q for s > with ds, d-^, and d^ all linearly independent 
at q. In particular, the phase function (1 — s — t)^ + Y ■ fi parametrizes (Aq, Ai) near the 
front face of Xq. 

Definition 13. We say that u G /'"(Xg; Ao, Ai; fi^ (Xq)) is a paired Lagrangian distribution 
of order m associated to Aq and Ai on Xq if it may be written as the restriction to Xq of 
some distribution in I"'{2X^; A^q, Al;n^2X^)). 

Definition 14. The class I^{Xq; Aq, Ai;^1^(Xq)) consists of the paired Lagrangian distri- 
butions that are supported away from the side faces of Xq. 

We similarly define the classes /™(Xq; Aq; fi^(Xo)) and /^(Xq; Ai; ^I^Xq)) as restrictions 
of distributions on 2Xq that are supported away from the side faces. 
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We typically shorten I^^iX^, Aq, Ai;n^X^)) and I^iX^, A;n^X^)) by suppressing the 
half-density factor. 

The following proposition follows immediately from Proposition 4.1 of |18] . 

Proposition 15. Ifu G J™(Xq; Aq, Ai), thenWF{u) C AqUAi. If B is a properly supported 
h-pseudodifferential operator of order on X^, then 

WF\B) n Ao = =^ Bue I^{X^; Ai \ ^Ai), and 

WF'{B) n Ai = =^ Bue /o™"'/'(X2; Aq \ 9Ai). 

Note that the orders on Aq and Ai differ by |. This can be seen by integrating by parts 
once in s. 

Just as in [18j, this gives us two symbol maps. For u G /™(Xq; Aq, Ai), 

a«(«)GC-+„/4(Al\5Al;^]^®Ll) 

a^'\u) G C^_i/2+„/4(Ao \ dA,; ® Lo), 

where the subscript indicates the degree of homogeneity and Li is the Maslov bundle over 
Aj. Admissible symbols are subject to a compatibility condition at dAi as in [18], but we do 
not need the explicit form of this condition here. We do, however, use the fact that if the 
principal symbol of u vanishes then u is one order better. 

For these classes to be well-defined, we must show that these classes are independent of 
the choice of extensions Aq and A^ (It was already shown in [12\ that the classes I^{Xq \ A) 
are independent of the choice of extension.) 

Lemma 16. The class /™(Xq; Aq, Ai) is independent of the choice of extensions Aq and Af. 

Proof. Suppose that Aq and Aq both extend Aq, while A^ and Af extend Ai. Let u G 
/q"^(2Xq2; A^, A^). We may find v G I^{2Xl, Ag, A^) with the same symbol on Aq U Ai. 

Because u and v have the same symbol there, u — v is order m — 1 on Aq U Ai. We now 
iteratively solve away principal symbols to find w with u — w of order — oo on Aq U Ai, i.e., 
such that u — w is a smooth function on Xq up to the front face, but supported away from 
the side faces, u — w is then clearly the restriction of an element of /q'°°(2Xq; A^, A^), proving 
the claim. □ 

The following proposition follows easily from Proposition 5.4 of |18j : 

Proposition 17. Suppose A is a properly supported b- differential operator of order m on 
Xq and is such that cr{A) vanishes on Ai. If u e oj^/q (Xq; Aq, Ai), then Au = x^h + x^g, 
with h G /q"'+''~'/'(X2; Aq) and g G I^+^^-^X^, Aq, Ai) . 

Proof. We start by noting that the action of Diff^(XQ) (extended to 2Xq) commutes with 
restriction to Xq. This follows from the observation that xdxX{x>o} = 0, where X{x>o} is the 
characteristic function of a half-plane. 

We start by writing u = x^ulx^, where u G /q (2Xq; Aq, A^). We apply Proposition 5.4 of 

[T8] to u, giving us Au = h + g. 

In particular, we have that Au = x^h + x^g + [A,xP]u. The operator A is a b-differential 
operator, so we know that [A, x^] = x^B, where i? is a b-differential operator of order m — 1. 
We thus have that Au = x^h + x'^{g + Bu), where h and g are the restrictions of h and g to 
X^. g,Bue Jo+'"~^(X^; Aq, Ai), proving the claim. □ 
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We also need another statement from Proposition 5.4 of [18], namely that the principal 
symbol of the distribution g above is given by the action of the Hamilton vector field of the 
principal symbol of A on the principal symbol of u. We only require this on the interior of 
the double space and so we do not prove it here. 

Definition 18. We define the normal operator of a paired Lagrangian distribution (or a 
Lagrangian distribution, as in [12]) to be the restriction of its Schwartz kernel to the front 
face ff . 

Note that this restriction is well-defined by wavefront considerations and the transversality 
of Aq and Ai to ff. Moreover, by considering the model form (fT2|) . we find that the restriction 
is a paired Lagrangian distribution in the class /"'(ff ; Aq, A^). 

Just as in [16], we have a short exact sequence. 

Proposition 19. The normal operator induces exact sequences 

(13) ^ xl^iX',; Ao) ^ I^iXl, Ao) ^ I^{S; A°) -> 
and 

(14) ^ xl^iX',; Ao, Ai) -> /-(X^; Ao, Ai) -> /™(ff ; A", A?) 

such that for any differential operator P G Diff™(X) and any paired Lagrangian distribution 
u, 

(15) Np{Pu) = Np{P)Np{u). 

Proof. The exactness of the sequences f|T3l) and f|T4|) follows from Taylor expansion in x near 
ff. 

The proof of f lT5|) is identical to the one in [16]. □ 

Also note that because LC is an embedded submanifold, elements of /o^^^^(Xq; Ai) may 
be identified with distributions of order k — 1/2 conormal to LC. 

9. Another blow-up 

In order to understand the solutions of the transport equations, we must introduce another 
blow-up. Because our solutions eventually have differing asymptotic behaviors along the light 
cone and on the interior of the light cone, we blow up the boundary of the light cone (i.e., 
at the left and right faces). This blow-up always makes sense locally near the front face, but 
to make sense of it globally we use assumption (A3). 

Because the boundary of the light cone meets the corner lf+ nrf_, the order in which we 
blow up the two submanifolds LC fl lf+ and LC nrf_ matters. We deal with this situation by 
blowing up the submanifold LC fl If + fl rf _ of the corner. We include this discussion primarily 
for completeness, as this piece of the construction is unnecessary as long as we restrict to 
forward-directed data. Indeed, the product of the pullback of a forward-directed function 
and any tempered distribution on Xq will vanish to all orders at this new face. 

We first define what we refer to as the intermediate double space. This is the space on 
which we solve the transport equations for the conormal singularity. If we restrict to data 
supported away from YL, this space suffices for our entire construction. 

Definition 20. We define the intermediate double space Xq^^ = [Xg, LC fl rf, LC fl If]. 
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This is a new manifold with corners. We will call lcf_|_ the lift of LC fllf and lcf_ the lift 
of LCnrf. 

Though we may think of this new manifold as an invariant way of introducing polar 
coordinates near LC fl If, projective coordinates are more convenient for our applications. 
Near lcf+ and ff+ but away from lf+, we may use coordinates p/s, s, and x. Similarly, near 
lcf_ and ff_ and away from rf_, we may use p/s, s, and x (and the remaining z and y 
variables), where s = x/x. 

Away from the front face and the corner lf+ fl rf_, this is just the blow-up of an intersection 
of two hypersurfaces {p = 0} and {x = 0} (or {x = 0}). Near lcf_ but away from lcf+ and 
fT_|_, p/x and x are valid coordinates, while near lcf+ away from lcf_ and fF_, p/x and x are 
valid coordinates. Near their intersection but away from lf_|_ and rf_, p/{xx), x, and x are 
vahd. 

Because lcf+ and lcf_ intersect, the order in which we performed the blow-up matters. 
Near the interiors of the new faces the two spaces are locally diffeomorphic, but they are not 
globally diffeomorphic. This is relevant when understanding the behavior of our parametrix 
near intersections of these faces, so we instead work on what we call the full double space. 
In this space we first blow up LCnlf+nrf_ and then perform the other two blow-ups. 
Performing this first blow-up separates the lifts of LCnlf+ and LCnrf_, so the order in 
which we perform the other two blow-ups is now irrelevant. 

Definition 21. The full double space is = [X^ LC n lf+ n rf_, LC n lf+, LC n rf_]. 

This is another new manifold with corners and is the space on which the final parametrix 
lives. We call the additional corner scf, for scattering face. We give it this name because 
this face is related to Vasy's construction of the scattering operator in [21]. Near the interior 
of scf nlcf+, we may use coordinates given by p/{xx), x/xt, and x, while near the interior 
of scfnlcf-, we may use p/{xx), x/x, and x. We again emphasize that this blow-up is 
unnecessary if we restrict to data supported away from y_. Figure H] depicts a neighborhood 
of ff+ in Xq (or, indeed, in Xg J. 

Ao 

y,y 

Figure 4. The double space Xq near ff+. 

Up to a smooth nonvanishing factor, the Jacobian determinant of the blow-down map 
Xqj — )• Xq is given by ricf^ricf_, so sections of the bundle Q2(^Xq) lift to sections of 
{rici+ricf^y^'^^^ {Xqj-) . Similarly, sections of ^1^{Xq) lift to sections of rscfrl^^^rl^^_Q,^ (Xq) . 

One can see that the fundamental solution does not lift nicely to Xq by considering 
the forward fundamental solution for the Klein-Gordon equation on the half-space MJl = 
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(0, oo)s X ^ equipped with the Lorentzian metric — ^J' ^ . Finding the forward funda- 
mental solution requires solving the equation 

P{X)u = 6{s-l)S{z), 
u = for s > 1. 

This solution can be constructed explicitly by taking the Fourier transform in the z variables, 
which transforms the equation into a Bessel ordinary differential equation for each value of 
the dual variable to z. By applying a stationary phase argument to the inverse Fourier 
transform of this family of Bessel functions, one can see that the asymptotic behavior of the 
forward fundamental solution is qualitatively different from its behavior on the interior of 
the light cone, justifying the blow-up. 

10. Outline of the construction 

We present now an outline of the parametrix construction, which consists of five steps. 
At each step we construct an approximation of the fundamental solution that captures the 
"worst" remaining singularity and yields an error term that is less singular. 

We simultaneously construct a distribution that is both a left and a right parametrix for 
our operator. In other words, we construct a distribution K such that P{X)K = I + Ri and 
P(A)^i^ = / + R2- Here the subscript R denotes the operator acting on the right variables 
and -P(A)* is the transpose of -P(A) as in equation ([5]). The distributions -Ri and R2 are 
smoothing operators that vanish to infinite order at relevant boundary hypersurfaces and 
should be considered "negligible". 

First we show that the construction in [18] is valid in a neighborhood of the diagonal, 
smooth down to the front face of Xq. The Schwartz kernel of the identity is a distribution 
conormal to the lift of the diagonal to Xq, and so the fundamental solution of -P(A) is a 
paired Lagrangian distribution associated to the lift of the diagonal and its fiowout via the 
Hamilton vector field, uniformly down to the front face. This is the most singular piece of 
the distribution but is the most standard portion of the construction. After cutting off this 
distribution outside a neighborhood of the lifted diagonal, this step yields a remainder term 
that is Lagrangian on the interior, uniformly down to the front face, and supported away 
from the side faces. Our assumption (A3) (which is unnecessary if we are only considering 
the construction in a neighborhood of y+) guarantees that the Lagrangian remainder term is 
in fact a conormal distribution associated to the projection of the fiowout of the light cone. 

In the second step of the construction, we solve away the conormal error term to obtain 
a remainder that is smooth on the interior of the double space. Because the light cone is 
characteristic for -P(A), this reduces to solving a sequence of transport ordinary differen- 
tial equations. These may be solved on the interior, but solving each successive transport 
equation on Xq causes a loss of an order of decay at the boundary. In other words, with 
each successive improvement of regularity at the interior, the behavior of the solution at the 
boundary becomes worse. 

We instead solve the transport equations on the intermediate double space Xq^. We solve 

these equations on Xq^ rather than the full double space Xq because transport equations 
have a slightly simpler form here. Indeed, the main modification needed to handle de Sitter 
space is to solve the equations on Xq rather than Xq^. This step then yields a remainder 
term that is smooth on the interior and poly homogeneous at the light cone faces lcf+ and 
lcf_ and the scattering face scf . 
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The third step solves away part of the polyhomogeneous remainder at the hght cone faces. 
We write down a formal series expansion for the solution and solve term by term, which 
yields a sequence of ordinary differential equations. The numbers s±(A) are the indicial 
roots of these equations near the side faces of Xq, so x^"^^^^ solve the equations to first order. 
In other words, the fundamental solution has an expansion at the side faces starting with 
3;«±{^) The remainder term from this step vanishes to all orders at the light cone face and 
has an expansion at the side faces. 

In the fourth step we solve away the remaining error term at the front face and side faces. 
This reduces to another formal power series calculation, and was carried out in [23]. The 
distribution we obtain by adding the distributions found in the four steps then solves the 
left (right) equation up to an error term that is smooth and vanishing to all orders at the 
"plus" ( "minus" ) boundary faces. 

In the fifth and final step, we remove this last error term via an argument showing that 
the fundamental solution itself must be in our class of distributions. 

11. The diagonal singularity 

The aim of this section is to solve away the diagonal singularity of the fundamental so- 
lution, leaving us with a Lagrangian error. The method here is similar to the one found in 
|12] . We use the normal operator to solve away the paired Lagrangian singularity at the 
front face up to 0{x°°), which allows us to invoke the construction in [T8] . 

Proposition 22. Given f G x^/q (Xq; Aq), there is a paired Lagrangian distribution u G 
x%"^/^+''(X2;Ao,Ai) such that 

p{X)u-fe~x^'l,"^^\xl-A,). 

In particular, if f is the Schwartz kernel of the identity operator, then 

f = 6is- l)5(^)i-"/V G Ao), 

and there is a paired Lagrangian distribution u G x~"'/^/o^^^^(Xq; Aq, Ai) such that P[\)u — 
fex-I^I^''\Xl-K{). 

Proof. On the interior of the manifold we may use the construction in [18], so we localize 
near the front face. Because we are only considering distributions supported away from the 
side faces for now, we do not need separate arguments for dealing with ff+ and ff_. 

We start by fixing a cutoff function x ^ C'°°(Xq) such that x = 1 in a neighborhood of 
the diagonal, but supported away from the side faces of X^. We also note now that -P(A) 
commutes with multiplication by x. 

Consider now the normal operator Np{P{X)) acting on the fiber ffp over the front face. In 
the projective coordinates (s, z,x,y), this is given by equation (|8|). 

It is easy to verify that Np{P{\)) satisfies the assumptions required to apply Proposition 
6.6 of [ig. We may thus find k G /^-^/^(fr; A[], A?) such that Np{P{X))K{p) - Npix'Pf) G 
C°°{S). Multiplying k, by the cutoff function x gives 

Np{PiX))ix<p)) - Np{x-^f) G /J-'/'(fr; A?). 
By appealing to the short exact sequence (1141) . we may find some paired Lagrangian 
distribution uq G 5;P/o~^''^(Xq; Aq, Ai) such that N{x~Puo) = x'^xi^- is characteristic 
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on Ai, SO we may appeal to Propositions [T71 and [T9] to write 

P{X)Uo - f = Vo + Wo, 

where vo E xP+^I^{X^; Aq) and wo € xPI^'^^^X^; A^). 

The distribution vq is supported away from the side face, so x~^Vo G x^Iq^Xq] Aq). We 
may now again use equation ( !T3|) . Proposition 6.6 of [18], and equation (fT4|) to find mi G 
£P/J"^/^(X2;Ao,Ai) such that 

N{P{X))N{S:-Pui) = N{xx-^vo). 

In particular, we may then again appeal to Propositions [T7] and ^\ to write 

P(A)(x(mo + XMi)) = V1+W0+ XWi, 

with t;i G xf+2/fc(x2; Ao) and G xPIo'^^'^iX^, Ai). 
Iterating this process gives us distributions Uj such that 

(iV \ JV 

i=o / j=o 

with Vat G xP+^+i/o^(X2; Ao) and G xPJo"^/^(X2; Ai). 

We may now asymptotically sum the x^uj to find u G x^Iq '^^^(Xq; Aq, Ai) such that 

PiX)ixu) -f = v + w, ve x^I^iXl, Ao), ^ G x%'-'/'(X2; Ai). 

We now extend v by zero to 2Xq, and use proposition 6.6 of [L8\ to find 

uelt''^\x',;Ao,A^) 

such that P{X)u + v E C°°{2X^). In particular, P(A)(xS) - t; = w G ^^/^"^/^(X^; Ai). The 
distribution + x£t then satisfies 

(16) PiX){xu + xu)-fe x^ll-"\xl- Ai). 

□ 

We should note here that solving the transport equations for the symbols on Ai (i.e., 
when invoking Proposition 6.6 of [18]) actually fixes the behavior of the solution in both 
components of Ai \ 0. In particular, we may arrange it so that the symbol on Ai has an 
expansion in decreasing powers of (77 + iO) (here rj is the fiber variable in Ai = N* LC). This 
allows us to guarantee that our parametrix u and the error term in flTBl) are supported on 
the interior of the light cone. This is not surprising because the exact solution must also 
be supported on the interior of the light cone due to the finite speed of propagation for the 
wave equation. 

The following lemma is also useful: 

Lemma 23. When f = kj is the Schwartz kernel of the identity operator, then the distri- 
bution u constructed in this section is also a left parametrix for P, i.e., Pj^u — f = x~"/^r, 
where r is smooth in a neighborhood of the diagonal. Here P^ denotes the transpose operator 
for P acting on the right factor. 
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Proof. This is an application of the symbol calculus for paired Lagrangian distributions. 

We start by observing that {Pl — P^)k,j = because the identity operator commutes with 
any operator. In other words, we must have that 



Now let V = {Pl- PA) 

It G X "^^/q (JCq ; Aq, Ai). We know that Plv is smooth in a 
neighborhood of the diagonal, down to the front face, because P^ and commute. We now 
invoke the symbol calculus: 

= (T (Plv)] Ao\dA, = (r{PL)cr{v)\Ao\dA,, 

and (j{Pl) 7^ on Aq \ 9Ai, so cr(f ) = on this set. This fixes an initial condition for (T(f ), 
i.e., a{v)\QA^ = 0. We now use the form of the transport equation 

= a{PLv)\A^ = {iChl + c)o-(t;)|Ai, 

with initial condition 0, to conclude that (j{v)\ai = 0. We may use this argument for any /c, 
so we must have that v E x~'^^'^Iq°^{Xq; Aq, Ai), proving the claim. □ 

12. The transport equation 

We now wish to solve away the error from Proposition[22j We call this error r and note that 
r G rg/Q~^^^(XQ; Ai). By solving a transport equation for some finite time and multiplying 
by a cutoff function, we may assume that this error is supported in a neighborhood of the 
side faces. Viewed as a conormal distribution near ff+, r may be written as 

(17) / e''^'^a{s,-,e,x,y,T])dr]. 



s 

where p is a defining function for the light cone and a is a classical symbol of order {k — ^) + 
In fact, as noted in Section [TTl we may assume that the symbol has an expansion in 
decreasing powers of [rj + iO). Near lcf_|_ nlcf_, we must replace ^ with 

We first fix a defining function p for the submanifold LC. Assumption (A3) guarantees 
that LC is an embedded submanifold of X^. Because A^* LC is characteristic for P{X) and 
dp spans A^* LC, we must have that g{dp, dp) = at LC, i.e. 

(18) g{dp,dp)=pb, 

where 6 is a smooth function. 

In coordinates near dX = Y, we may write 

P(A) = {sd,r - (n - l)sd, + ^^^sd, + s'A, - A. 

Our ansatz is that u is polyhomogeneous at lcf+ and Icf _ and conormal to LC, and so we 
seek an expression of the form fll7p with a a classical symbol of order A; — | + = + | — 2 
with an expansion in powers of [r] + iO), i.e., of the form 

(19) « ~ 5Z / ^'^"(^ + ^0)''+^"'~^a,- d7]. 

This is because our error r from Section [11] is a Lagrangian distribution of order k — ^ 
associated to the conormal bundle of LC. LC is characteristic for P(A), so we expect the 
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solution of P{X)u = r to be conormal of one order better, i.e., of order A; — |. Moreover, 
at each step, we multiply the symbols aj by a compactly supported smooth function in | 
that is identically 1 near 0, which makes the singularity of (t] + iO)''~^^~'^~^ at superfluous. 
Note that we could equivalently insist on expressing our ansatz in powers of (f If = 0, 
then the top power seen here would be — 1 + 2 — ^ = 1 — f, the same powers seen in the 
construction of the fundamental solution of the wave equation on Minkowski space (in our 
convention, n is the total dimension of the spacetime). 

Lemma 24. Suppose u is of the form fW\) . and aj = djU, where u is a fixed nonvanishing 
section of Vl^ (-^o *) ■ ''^^ write 7j = /c + | — 2 — j, then near lcf+ away from lcf+ fl Icf _ we 
have 



oo „ 

P(A)m~^ / e'-s~^{r] + iQf^^^{-2ig{dp,sdaj)-i{dsp){n-?,-2^j + 0{s))dj) dr]u 
+ ^'^'(^ + '0)"' (^(^) + + 0{s))sds - (7, + - 3 - 27,) 

3=0 ^ 

(20) +]^(^-^ +0{s)sd, + 0{s)^~ajdr]u 



where 0{s) is taken to mean an element of sC°° and (s = |,z = are projective coordi- 
nates near the front face. 

Proof. We first show the result near ff_|_ nlf+. 

Write Uj = uju, where z/ is a fixed trivialization of Q2{Xq^). Say, for concreteness, that 
u = r^J^'^x^^^'^v. We then have that 



P{\){u,u) = (P(A)S,) V + rl/lx-/' ([PW,r;Ji'S:~-/']) n, 



We note also that near lcf+ but away from lf+, ricf_^ = s(l + as), where a is smooth, and 
so we may easily calculate this commutator: 

irlIlx-/')[P{X),r;Jl'x--/'] = (-1 + 0{s)) sd^ + 0{s)sd. + ^ " ^ + ^(^)) ' 



We now use this calculation to drop the density factor. We apply P{\) to our ansatz and 
use equation ( !T8|) . Integration by parts allows us to exchange powers of ^ for decreasing 



24 DEAN BASKIN 

powers of {r] + zO), as ^e'^^'^ = \dr^e^^^ . If we write 7_,- = + ^ — 2 — j, this yields 
P{\)u = (^P(A) + {-l + 0{s))sds + ]^[n-]^ + 0{s)^ ^ u 

oo „ 

~ 5Z / ^'^"(^ + ^^V'^^9{dp, dp)aj dr] 
+ £ / e'^^^iv + iOp+^f -2ig{dp,sd&j) 

— i{dsp){n — 1 — 27j — 3)5^ — i{dsp)dj + 0{s)dj^ drj 
+ / + ^0)^' + 2(7, + - (7, + l){n-l- 27, - 1) 

+ (7j + 1) - + ^ (^n - + 0(s)s9, + 0(s)s9, + 0{s)^~aj dr], 

where 0{s) is taken to mean an element of sC°°. We now use that g{dp,dp) = pb and 
integrate by parts to prove the first part of the claim. 

We finish the proof with a similar calculation, where we change our operator to and 
our ansatz to be one of the following two forms: 



oo 

U ~ 

j=0 

oo 

U ~ 

j=0 



oo „ 

J2 / e'^'^aj{r] + tOp dr], 

3=0 
oo „ 

^ / e'^''aj{r] + iOy^ dr], 



where s = |, and P|. is given by equation (|5]). We use here that ricf_ = s(l + as) in the first 
case and that ricf_ = x(l + ax) and ricf_^ = x{l + ax) in the second case. □ 

Lemma 25. Similarly, we may compute the right operator near Icf _ away from the corner: 

oo „ 

P(A)*jM = XI / ^''''(^ + ^0)^'^^ {-2ig{dRp, Sdndj) + i{dsp){n + 3 + 27, + 0{s))dj) dr] v 
j=o 

00 ^ 
+ ^e^f.(^ + i^rnPiXYn + (7, + + (7, + + 3 + 27,) 



1 

2 V" 2 



+ \{n-\] + 0{s)sd-^ + 0(s) ) Oj dr] v 
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We may also compute the behavior of the left and right operators near lcf+nlcf_; 
P{X)ur^Yl / e'^'^iri + iOy^^^- {-2ig{dp, xdttj) - i{d^p){n - 3 - + 0{x))aj) dr/u 

+ I e'^"!^ + ^0)^^ (P(A) + (7,- + 0{x))xd, - (7,- + l)(n - 3 - 27,) 
i=o ^ 



\\^^\]+ 0{x)xdy + 0{x) ) aj dr] u 



2 V 2 

00 



r 1 

^(A)flii = 5] / ^'^"('^ + iO)^^+'- {-2ig{dRp, xdRUj) + i{dip){n + 3 + 27,- + 0(x))a,) ^77 1/ 
+ e'^'iv + iOr^ ( ^(A)*^ + (7,- + 0{x))xd, + (7,- + l)(n + 3 + 27,) 

+ - ^ 2^ 0{x)xdy + 0(5;)^ ajdrju 

Because we wish to solve P{\)u = r up to smooth terms, we wish to iterativcly solve away 
the terms in the above expansions. The first transport equation we must solve is then 

-2g{dp, sddo) - {dsp){n - 3 - 270 + O(s))ao = Tq, 

where tq is compactly supported and comes from the inhomogcncous term. 
Because dsP is nonzero, wc may divide by it to obtain the transport equation 

+ sdsao - (77 - - 7o + O{s))ao - 



dsP '"^22'" ^ " 2d,p 

Note that 70 = /c + | — 2, so the coefficient of is just k — ^ + 0{s). 

Near the face lcf+, given by s = 0, we may use a parameter t along the light cone LC. 
The parameter t is then equivalent to s, so we may change coordinates to 

1 

tdftto + {k- -)ao = 0{t)ao + f. 

We note that the solution oq of this equation must have a polyhomogeneous expansion in t. 
In fact, we prove a more precise version of this statement. 

Lemma 26. Suppose that v solves the differential equation 

tdtv — {j — k + = t ■ c{t)v + b, 
where c is smooth in t, and b is polyhomogeneous in t with index set 

£j-i= !^(^-k+^ + l,z^ ■.leNo,i<lifl<j-l,i^j-lifl>j 

when j and b is supported away from t = for j = 0. Then v has a polyhomogeneous 
expansion in t with index set 

i(-k +^ + l,tY. I e No,i<Uf I <j,i^j if l> J 
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Remark 27. We may prove a similar lemma for solutions of the transport equations on the 
right, with appropriate modifications for the index sets. 

Proof. This lemma follows as a simple exercise in the b-cal cuius of Melrose (cf. [7] or [21]) or 
as an exercise in the theory of hyperbolic Fuchsian operators (cf. [2j,[13j, or [22]). Because 
the proof of this lemma is elementary, we include it here. 
We show this by constructing a formal power series solution. 

We start with the case j = so that near t = 0, b = 0. We seek a formal power series 
solution of the form 

oo 

Indeed, with this ansatz, the equation becomes 



-1) 

1=0 1=0 

where depends only on c(t) and the first /—I coefficients Vi (so g_i = 0). The coefficient Vq 
is fixed by the initial condition of the differential equation, and then the remaining coefficients 
may be found iteratively. 

We then sum this series with Borel summation to find a function with this power series at 
t = and the difference between this function and the solution of the differential equation 
vanishes to all orders at t = 0. 

For general j, we write 

j-l min{l,j-l) 

^ = E E but\\ogty 

1=0 i=0 

with a similar expression for v. A similar calculation then reduces the equation to 

oo m'm(l,j) 

1=0 i=0 

oo I 

1=1 i=0 

where qi-i^i depends on the function c(t) and the coefficients vi>i> where V <l — 1 oi V = I 
and i' > i. In particular, we may again iteratively solve for each coefficient and then Borel 
sum the result. □ 

We may now apply Lemma [26] to the first transport equation 

tdtdo - U - k + ^)do + tc(t)do = b(t), 

where c{t) is smooth in t. We find that is polyhomogeneous in t with index set Sq. 
Changing coordinates back to s tells us that oq is polyhomogeneous in s with index set Sq. 

Letting Qj be the operator acting on dj in the coefficient of {t] + iO)'^°~^ in equation (120|) . 
the jth transport equation is then 



+ sdsaj — {j — k -\ — )aj + 0{s)aj = —Qjaj^ 



dsp ' ' 2 
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By applying Lemma |26] again, we may conclude that aj is polyhomogeneous in s with index 
set Sj. 

Now let S = UjSj. By repeating the process above, we obtain conormal distributions Uj 
in /^^'^^^^"' (Xq LC) with symbols aj such that 



We now wish to asymptotically sum this expression to find u such that P{X)u — / is 



smooth on the interior of Xq^. 



Lemma 28. There is a distribution u E ^''^(Xq j; LC), supported on the interior of the 
light cone, such that 

oo 

i=o 

Proof. Each symbol {T]+iO)'^°~^aj is analytic (in rj) and exponentially decreasing in the upper 
half plane '^rj > 0, so the Paley- Wiener theorem tells us that each Uj is supported in the 
region > 0}. A standard Borel summation argument completes the proof. □ 

Putting our factors of rg from Section [11] back in, we have thus proved: 
Proposition 29. Given r G r^I^ ^^^(Xq;LC), there is a conormal distribution, 



u 



ph-g 

such that 

P(A)n-rG«^^(Xj,;lcf+), 
where u is supported away from the corner ff+ fl lf+ and T\ is given by 

•^i = |(j-A; + i,/):j,/GNo,/<J 

We may further arrange that the distribution is supported in the interior of the light cone. 

We observe here that an argument similar to the one given in the previous section shows 
that P^M — r is again smooth on the interior when we are constructing a parametrix for the 
fundamental solution. In particular, the symbols of u satisfy the transport equations for the 
actions on the right factor, and so the following proposition holds: 

Proposition 30. There is a conormal distribution u G rg"^^^-^/^'^/^(XQ LC) such that 

P{X)u-rer-/'Alt,{Xl), 
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Here the index families T and £ are given by 



Icf 4 



*lcf_ 



lcf+,L 



lcf_,L 



E, 



M+,R 



E, 



lcf„,R 



(j + 2'0:J,^eNo,/<j 



(J + 
(J + 



- n, /) : J, / e No, / < J 

-n-l,0:j,/GMo,/<j 
-l,/):j,/GNo,/<J 
-n,l): J, / G No, / < J 



Note that the decrease of —1 for two of the error terms in the above proposition come 
from the factors of - in the expression given in Lemma [23 The —n on the "minus" faces 
comes from taking the transpose of our operator. This is because we are using sections of 
the standard half-density bundle rather than the 0-half-density bundle. 

13. The light cone face 



We now wish to solve away the error from Proposition |30l We aim to solve away the left 
error at lcf+ and the right error at lcf_. We show only the left calculation here and observe 
that the right calculation is nearly identical. 

We again call the error r, which is in r^A^^g{XQ^), i.e., near ff+, r has an expansion of 
the form 

oo j 
j 1=0 

where Jj = j — k — i. We drop the power of x for now because -P(A) commutes with x. 
We first claim that this error lifts to be polyhomogeneous on the full double space Xq. 

Lemma 31. Suppose that r is a polyhomogeneous distribution on Xq^ supported in a small 
neighborhood of LC with index family J-". Then r lifts to a polyhomogeneous distribution on 
Xq with index family Q, where Q is given by 

Clef. = Fic{_, 

Cscf = -Flcf+ + -^lcf_ + 1- 

Remark 32. Here the notation F + 1 is shorthand for 

(a, /) G F + 1 if and only if {a -1,1) e F. 

Proof. The result follows because the two possible orders of the blow-up are locally diffeo- 
morphic near LC nlcf+ nlcf_. The extra 1 in Ggcf is because sections of VL'2{XQ^) lift to 
sections of Tscff^ 2 (Xq). □ 
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We may thus consider r as a polyhomogeneous function on Xq. 

We now proceed in two steps. The first is to solve the away near Icf + and the second is to 
show that it has the desired form at lf_|_. This only away from the scattering face scf , though 
the computation near scf is nearly identical. In this section there are many terms that come 
from differentiating our ansatz. We attempt to indicate the origin of the important terms. 

Note that the statement about the support in Proposition [501 means that rji is supported 
on the interior of the light cone and vanishes to infinite order at LC. 

Because we are working near lcf+, we first use projective coordinates {s,w = -,0), where 
p is a defining function for LC (as above) and 6 are the remaining variables. In these 
coordinates, derivatives of the function p appear as coefficients of dw in our operator P(A). 
In particular, sds lifts to sdg + {dsp)dw — wdw 

We again expect a polyhomogeneous expansion in s. In other words, we expect an expan- 
sion of the form 



(21) «~^5^s^^(logs)V 

j=0 1=0 

where Uji is regarded as a function of w, 9, x, and y. 

Write Uji = Ujiv. We again note that we may again take ri^f^ equivalent to s in this region, 
giving us an extra —\sds — |(n — i) in our operator. Applying P(A) to our ansatz yields 



P{\)u r-J2Yl '''' [-9idp, dp) - 2w{dsp) + w^) {dlu^i) 

j 1=0 

+ 2'jj){dsp) - (n + i)w^ d^^Uji 
(n - 1 - 7j)7j + ^{n - ^) - ^ji 



where 



X.ji = 2{l + l){dsp - w)d.^ -{1 + l){n - ^ - 27,), Yi = {l + + 2). 



The constants above come from the operator and from sds landing on the powers of s. The 

'z,p){dz,p)dl 



leading 3^ terms come from the {dsp)dw terms when we lift sds and the h^\dz^p){dziP)d'f 



term in the lift of the Laplacian in z. 
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Note that because LC is characteristic for □, we must have that g{dp, dp) = ph. p 
so by replacing -Bj+i with = -Bj+i — wbdl^, we may write 

j 



1=0 



where 



Q, = i-{d,p) + w) wdl - {id,p){n + i - 27,) - (n + j d, 

-{n-l- 7j)7j- - ^(n -\) - ^■ 

We thus wish to solve a sequence of inhomogeneous ordinary differential equations QjUji = 
Vji, where 

where all terms are supported in {w > 0}, vanishing to infinite order at w = 0. (We know 
already that rji has this property, and we show at each step that Uji does as well.) 

It is clear that is a regular singular point of the differential operator Qj, so the solutions 
of QjV = have formal power series expansions at w = with first term given by w'^' . Here 
Pi are the roots of the indicial equation —p{p — 1) — (n + | — 2'jj)p = (see, for example, 
P), i.e., 

Pi =0, P2 = 2j —2k — n 

Here we have used that 7^- = j — A; + |. Standard ODE techniques (i.e., variation of param- 
eters) then give us solutions Uji to QjUji = rji in terms of a basis of solutions for QjV = 0. 
Moreover, because rji vanishes to all orders at w = 0, we may also guarantee that Uji is 
supported in {w > 0}, vanishing to all orders at w = 0. Indeed, if vi and V2 are a basis for 
the solutions of QjV = 0, then 

Uji{w) = -vi{w) / — -xr^dw +V2{w) / — '—-—dw. 

Because rji vanishes to all orders at w = 0, (and Vi are bounded by w^^ for some A^), the 
integrals make sense and vanish to all orders at w = 0. Multiplication by Vi then preserves 
this property. 

We may thus solve these equations and now wish to consider their asymptotics near the 
corner lcf+nlf+. Near this corner, the coordinates {s,w,6) are invalid and so we must use 
the other set of projective coordinates (p = sw, W = w~^, 0). In these coordinates, sdg lifts 
to {dsp)dp + Wdw — {dsp)W'^dw- The top order terms below then come from the Wdw and 
the constants come only from the operator. 

Because polyhomogeneous distributions are independent of our choice of coordinate sys- 
tems, we may also express our ansatz fl2T|) in terms of these coordinates {p,W,6). In this 
case, expansions in s are equivalent to expansions in p, so our ansatz here has the form 



5^5^p^^(l0gp)S7. 



U 

j 1=0 
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In the computation that follows, the important point is that S-i-(A) are the indicial roots of our 
operator -P(A), and that this behavior is dominant away from LC because the fundamental 
solution is a smooth solution of the homogeneous equation here. 
Apply P(A) to such an ansatz to see 

P(A)«- J2 f^{^ogpy[{l-q2W)W'd^wV,i 

- (n - 2 + qiW)WdwVji - (A + qoW)vji 

where we have (via a similar, but more tedious, calculation) 

q2 = 2{dsp)-{dspfW+\d,p\lW, 
q, = -{n + 27, - i){dsp) + (27, - 2)g{dp, dp)W, 
qQ = {n- 2)^j{dsp) - - l)g{dp, dp)W, 

while 

Aj_,,i = 2h^\dMid.,p)W'de^dw - 2^,h^''{d,^p){dMW\ 

-iDgp)w^dw-iAn-,p)w, 

= (/ + !)[ {-2{d,p) + 2g{dp, dp)) Wdw 

- {dsp){n - 2) - (27,- - 2)g{dp, dp)W] , 

and 

Bj_2,i = {/^hep)Wde, - }^\d,,ep){dMWde,de, 

= -2{i + i)h}\d,,p){dMwde^ + (z + i){n,p)w 

S^^2,/+2 = -{I + + mdp, dp)W. 

We must still take into account the density factor here. In this region, ricf_,_ is equivalent to 
p. A similar computation to the one above shows that this leads to two types of terms. One 
of these is 0{W)Wdw^ while the other is 0{W). This means that they may be absorbed 
into qi and go- Let us now write 

Qj = (1 - q2W)W^dl, - (n - 2 + qiW)Wdw - (A + q^^V). 

We wish to solve a sequence of transport equations (which are really the same equations as 
above written in these coordinates) given by 

Qj'fjjl = Tjl - WAj_i^iVj_i^i - ,+i^j„i,;+i - WBj_2,lVj-2,l 

- WB'j_^ i^-^Vj^2,l+l - W^5j'-2,«+2^i-2,i+2 = f'ji. 

The solutions of these transport equations are polyhomogcncous conormal functions. Here 
the indicial roots of Qj are s±(A) — ± y^ *-""^'' + A. More precisely, we have 

Lemma 33. Suppose that Uji solves QjUji — f'^i as above. Suppose first that s+(A) — s_(A) ^ 
Z. Suppose that f'-^ is polyhomogeneous conormal at lf+ with index set F + 1 — (F+ + 1) U 
{F~ + 1), where 

i^^ = {(s±(A)+m,0) :meNo,}, 
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then Ujl is polyhomogeneous with index set F = U F~ . 

//s+(A) — s_(A) = N E "Z, then if f'-^ is polyhomogeneous with index set F + 1, where 

F = {{s_{\) + mj) : m eNo,l = for I < N,l = 1 for I > N}, 

then Ujl must be polyhomogeneous with index set T . 

Remark 34. Note that the index set above means that if s+(A) — s_(A) ^ Z, then Uji G 

s+(A)-s_(A)-l 

u^i - ^ ly^- W+^mJ.^) g (1 + log iy)C°°. 

m=0 

Proof. We may again construct a formal power series solution (or apply Lemma 5.44 of 
Melrose [21j). We omit this here because it has been described in detail already. The key 
point is that the operators -Bj-2,i, -Bj_2,i+i! and -Bj_2_;+2 are all elements of 

VTDiff^ (products of W and differential operators tangent to the boundary). This ensures 
that the terms we are solving away vanish to one order better at ly = 0. □ 

Asymptotically summing YlYl P'^^ i^'^S pY'^ji then solves away the error r at lcf+. Note 
that because s±(A) are the indicial roots of N{P{X)), we in fact have that the error term 
vanishes to one order better. 

When we are constructing a parametrix for the fundamental solution, we may perform the 
same construction on the right and the left. We may also add a distribution solving away 
the right error at lcf_. 

Remembering our factors of rg now, we have now proved the following proposition. 

Proposition 35. Given ri G rfj"^^^p^^(XQ) and r2 G rg"^^^p^(XQ) as above, we may find 

a smooth function u G '"^"''^.^^^(Xq) vanishing outside the light cone such that P{\)u — ri 
vanishes to all orders at lcf+ and is polyhomogeneous with index family Qi, while P(A)^m — r2 
vanishes to all orders at lcf_ and is polyhomogeneous with index set Q2- Here we have that 

(22) Fif^ = Gif+,2 = {(s±(A) + m, 0) : m G No}, Fi^^ = Ficf+,L, 

Fri+ = Grf+.i = {i-n + s±(A) + m, 0) : m G Nq}, Fm_ = EM_^n, 

G'icf_,i = -Ficf_ — 1, ^if+,1 = -^if+ + 1) 

Clef +,2 = -^lcf+ — 1, Crf+,2 = -^rf+ + 1, 

Clcf+.l = G'icf_,2 = -^lf+,L = -f'rf-.L = -^If+.-R = -^rf+,i?, = 0) 

Fscf = Escf,L = Esci,R = Gscf,i = G'scf,2 = {{—n + j, • ^ ^ i}) 
z/s+(A) — s_(A) ^ Z. If s+ — s_(A) = X G Z, then the index sets become 

Fif^ = {(s_(A) + mj) : m ^^0,1 = for m < N,l = 1 for m> N}, 

Frf^ = {{-n + s_(A) + m, /): m G No, / = /or m < X, / = 1 for m > N}, 

with corresponding changes for Gu, and Grf+ • 
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14. The front face 

We now wish to solve away the errors on the front face, which we again call r, from the 
previous step. We show how to solve away the error term at ff_|_ for the operator acting 
on the left, and the corresponding calculation at ff_ for the operator acting on the right 
is nearly identical. We now suppose that r is the error term from Proposition for the 
operator acting on the left. 

Because r vanishes to all orders at lcf+, we may blow down lcf+ to solve away r. In this 
view, x~^r is smooth on ff+, supported inside the light cone, and has an expansion at lf+ of 
the form 

r ~ s"^ (logs)Vj7/i. 

The generic case here is that / = when s+(A) — s_(A) ^ Z. We again drop the powers of x 
from our notation. 

We wish to solve this error away with a function of the same form: 

u~ ^ s"^(logs)'Mj7/^. 

Applying iV(P(A)) to such an ansatz yields 
(23) 

NiP{\))u ^"^(log^)' [«'%^ + 2%(/ + + (/ + !)(/ + 2)u,j+2 

j I 

-{n - l)ajUji - {n-l){l + + ^zUj-2,1 - >^Uji] ^ 

= J2Y1 sY [- ((n - 1 - a,)aj + A) uji - (n - 1 - 2a,){l + l)uj,i+i 

j I 

+ {l + !)(/ + 2)M,y+2 + A,Uj^2,l] /i. 

The coefficient of Uji in this expression vanishes precisely when aj = s-i-(A), so we may solve 
away the error to all orders at lf+ because the expansions of r begin at s-i-(A) + 1. 

Proposition 36. For s+(A) — S-(A) ^ Z, there is a smooth function u on ff+, polyho- 
mogeneous at lf+ with index set F^^, where F\{^ is defined in equation [2^) . such that 
N{P{X))u — x~^r vanishes to all orders at lf+nff+. An identical statement holds when 
s+(A) — s_(A) G Z, with F\i^ replaced by F\i^. 

Proof. This is clear from the expression fl23l) for N{P{X))u. □ 

On the front face, we are then left with an error term tq vanishing to all orders at lf-|. fl ff + 
and supported inside the light cone. Because ff+ is an asymptotically de Sitter space, we 
may now use Corollary 3.6 of Vasy ([2^) to find a smooth function on ff+, vanishing to all 
orders at lf+, such that N{P{X))u = x'^tq. 

By iterating this construction (and extending it to the interior of lf+ as well), we may find 
u G a^^"4^g(^o ' ^^^^ that P{X)u — r vanishes to all orders at ff+ and lf+. By gluing two 
of these functions together, we may simultaneously solve the left error near the "plus" faces 
and the right error near the "minus" faces. This proves the following proposition. 
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Proposition 37. We may solve away the left error at ff+ and lf+, and the right error at 
rf_ and ff_ with a function u G ^^^(Xq), where H is given by 

Hm+ = Hici_ = 0, = for the other index sets. 
The remaining error terms are in •A'^hg'^ {Xq) , where JCl and JCr are given by 

K\i+,L = -f^rf_,iJ = -^lcf+ = Kicl_ = 0, ^\i+,R = -^lf+' 

15. The full parametrix 

We now take the various pieces of ttie parametrix constructed in Sections [TT], [T21 [13 and 
[T^to construct a parametrix for the fundamental solution of P. 

Putting together the results of Propositions [221 [30l [351 and[371 we have proved the following 
theorem: 

Theorem 38. Suppose that X is an asymptotically de Sitter space, satisfying assumptions 
(Al), (A2), and (A3). We may find a left parametrix K such that P{\)K = I + Ri and 
KP{X) = I + R2, where the Schwartz kernels of Ri and R2 are smooth on the interior of 
Xq and are polyhomogeneous with index families £l and £r on X^. We may write K = 
Ki + K2 + K3, where Ki is supported in a neighborhood of the diagonal, K2 is supported in 
a small neighborhood of the light cone LC away from the diagonal, and all three pieces are 
supported on the interior of the light cone. Moreover, 

(24) Gr-"/ V^'(X2;Ao,Ai), 

K, e r-"/Ui(Xo^). 
The index families £l, and T are given by 

Fm+ = {{j + ^J) -.1 <J,J e No}, Fm_ = Ficf+ - n, 

Fu+ = ^if+,ij = {(s±(A) + m, 0) : m G No}, ^rf_,L = i^rf_ = i^if+ - n, 

Em_,L = -^Icf- — 1, -Elcf+.iJ = -^lcf+ — 1, 

Fscf = Escix = Esci,R = {{-n + : I < j}, 
Es.,L = Es^^R = I (^-^ + m, 0) : m G No} , 

Elci^,L = -Elcf„,/? = -Slf+,L = E^{_ji = -Efr+,L = Eg_ji = 0. 

//s+(A) — s_(A) G Z, then we must modify F^^ and Frf^ (and the index sets depending on 
them) as described earlier. 

Remark 39. As observed earlier, this theorem holds without the assumption (A3) as long as 
we are willing to multiply our distribution by a cutoff function supported in a neighborhood of 
the front face. In this case, the extra blow-up to obtain the scattering face scf is unnecessary. 

Because our remainder terms lose one order of decay at the light cone faces, we lose an 
order of decay there when we pass to the exact fundamental solution. The following is a 
precise statement of the main result (Theorem [1]). 
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Theorem 40. The exact forward fundamental solution is in this class of distributions, 
but with index sets F{^^^ = Ficf^ — 1 and F{^^_ = F\ci_ — 1. In other words, we may write 
= Ki + K2 + with 

(25) K,er~^^'l~'/\x',;Ao,Ai), 

Here the index family T is given by 

Ficu = {{j -^,l):l< J, 3 e No}, = Fj^^ - n, 

Fu^ = {(s±(A) +m,0):me Nq}, Frf_ = F^^ - n, 

F,,f = {{-n-l+j,l):l<j}, 

with modifications to F^^ and Frf_ when s+(A) — s_(A) is an integer. 

Remark 41. If we instead adopt the convention that the 0-densities are flat, that 

fix, y)6{x - x)6{y - y) dg{x, y) = /(x, y), 

and we write i^T as a section of the pullback bundle of x X), then the index sets 

change somewhat. Indeed, if = Ku, where z/ is a nonvanishing section of the pullback of 
^Q^{X X X), then the same theorem holds, but with index sets 



-^lcf+ 


— -^Icf- 


_ = {{j-l,l):l<j,jeNo}, 


-^lf+ 


= Frf_ 


= {(s±(A) +m,0) : m E No}, 


Fscf 


= {(j - 


-l,/):/<j,jGNo}, 




= Fs_ 


= {(j,0):jGNo}. 



Proof of theorem. E^ is the forward fundamental solution, so if / is a compactly supported 
smooth function on X, then P{X)E+f = f and E+P{X)f = f. Moreover, continuity allows us 
to extend this to any forward-oriented distribution. In particular, if / is any smooth function 
on the interior of X vanishing to all orders at that is also a tempered distribution on X, 
then 

P{X)E+f = f 
E^P{X)f = f. 

Let K be the parametrix for £"+ constructed in Theorem [381 If / is a smooth function on 
X, vanishing to all orders at F , then Kf vanishes to all orders at Y_ because K is identically 
zero in a neighborhood of lf_ and the lift of / vanishes to all orders at lcf_, ff_, scf, and 
rf_. A similar argument applies to Rif and R2f ■ We may then write 

[Kf) = E^P{X)Kf = E^f + E+RJ 
(Kf) = KP{X)E+f = E+f + R^E^f. 

In particular, 

E+ = K- KRi + R2E+R1. 
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We then observe that the error terms KRi and i?2-S+-Ri have the desired properties, finishing 
the proof. □ 

16. Modifications for de Sitter space 

De Sitter space does not satisfy assumption (A3) because the projection of the forward 
fiowout of the hght cone from a point at past infinity intersects itself at future infinity 
(though not on the interior of the spacetime). We briefiy discuss the modifications to our 
construction needed for de Sitter space. The most important modification is to solve the 
transport equations on Xq rather than Xq^. (Indeed, we could have done this from the outset, 
but chose not to because solutions of the transport equations are easier to understand on 

Observe that the construction detailed above works without difficulty away from the corner 
If 4. nrf_. As mentioned earlier, LC intersects itself in this corner. This intersection is given 
by {(0, y, 0, —y) : y G Y}. Near this intersection, we may write p = x + x — \y + y\, plus terms 
vanishing to higher order at the boundary. Note that the function p is no longer smooth at 
this intersection because ||/ + ^| =0 there. 

We resolve this singularity by blowing up the submanifold LCnlf+nrf_ as in Defini- 
tion [211 Although this submanifold had codimension 3 when X satisfied (A3), it has codi- 
mension n + 1 here. Indeed, this blow-up is almost the same as the one in Section |4] that 
defined Xq. The function p now lifts to be smooth on this new space [Xq, LC fl lf+ fl rf_]. 
Indeed, after the blow-up we may write 

(26) p= {1 + s) -\z\+0{x) 

near scf fl If^. This is now smooth near p = 0, and so we may blow up LC fl lf+ and LC fl rf_. 
to obtain Xq. 

We must also modify the manner in which we solve the transport equations. The behavior 
at lcf+ and lcf_ may be obtained in the same way as in Section [121 but the behavior at 
scf requires a slightly different approach. In Section [9l we showed that the symbol of the 
conormal distribution was polyhomogeneous at scf by constructing it on Xq ^ and lifting it 

to Xq. Because p is not smooth on Xg when X is the de Sitter space, we cannot solve the 
transport equations on the intermediate double space Xq ^ up to the corner and instead we 
must solve the transport equation along lcf_ and lcf+. 

Solving the transport equation on these faces requires using the semi-explicit form of p. 
The terms where the operator lands entirely on p in equation f[20|) can no longer be ignored. 
Using the form f[26l) for p, we observe that these now contribute a constant term to the 
equation. Because Dp = s — {n — l)s + s'^{n — 1) + 0{x) and dsP = 1 + 0(5;), 

A computation in the same spirit as those in Section [T2] shows that the symbol of K2 is 
also polyhomogeneous at scf with index set 

Fsci = {{-n + m,0) -.m e No} . 

Note that this is the same index set we found before. 
The rest of the construction proceeds without change. 

The modifications to the construction for de Sitter space correspond to allowing the lo- 
cation of the pole p in P{\)u = Sp tend to past infinity. If we require that the point p is 
uniformly bounded away from past infinity, no modification is necessary. 
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In [23] and [25], the authors do not consider sending this pole to past infinity. Our 
unmodified construction recovers a shghtly weaker version of the results of these authors 
when the pole is in the interior but the modified construction also describes the behavior of 
the fundamental solution when the pole is at past infinity. 

17. POLYHOMOGENEITY 

The aim of this section is to prove Theorem [2], which was stated in Section [1] 
We begin by considering the maps /3l and where (3^^^ are given by composing the 
blow-down maps with projections onto each factor, as in the diagram here: 




X X 
We require the following four lemmata. 

Lemma 42. The maps o-nd Pi\i^Qint are b-fibrations . 

Lemma 43. The fibers of and /3l : LC™* — )■ X are transverse to LC. 

Lemma 44. For V a b-vector field on X , there is a b-vector field V on Xq such that 
= V and V is tangent to LC. 

Lemma 45. Suppose that M and N are manifolds with corners, F : M ^ N is a b-fibration, 
H is a boundary hypersurface of M, and F{H) = N Suppose also that K G I^{M;H) is 
polyhomogeneous at the other boundary hypersurfaces of M (satisfying the hypotheses of 
Lemma\9\j. Then F^K is a polyhomogeneous distribution on N . 

Proof of Lemma Here LC™* is taken as an open manifold away from the cone edge near 
the diagonal. K2 is supported away from the cone point, so we may restrict our attention to 
this region. 

That and /J^^Ilc™' b-fibrations (defined in Section [5]) follows from a more general 
statement: ii F : M ^ N is a. b-fibration, and Z C DM is a p-submanifold, then F o p : 
[M] Z] — )■ iV is a b-fibration. In particular, the blow-down map is a composition of b- 
fibrations. □ 



Proof of Lemma 43 Because Ai is the fiowout Lagrangian in the right factor, its intersection 
with the fibers xs = Xq, y + xz = is the fiowout of the light cone with cone point at (xq, yo). 
In particular, this is an embedded submanifold of X, and so the intersection is transverse. □ 



Proof of Lemma 44 This is really a consequence of Lemma Indeed, we may choose a 
local basis of vector fields given by di, and Vj, where di, is tangent to the fibers of (3^ and Vj 
are tangent to LC. Then = because di, is tangent to the fibers of /3l, and so we 

may choose a lift of the vector fields so that the diy component vanishes at LC. □ 
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Proof of Lemma \J5\ We choose K"^ polyhomogeneous on M, supported away from H such 
that K"^ are uniformly bounded in I^{M;H) and converge to K in I^{M;H) for any 6 > 
0. The pushforward theorem of Mehose (Lemma [9] of this paper) tells us that F^.K'^ are 
polyhomogeneous with fixed index set. 

We claim now that the F^K'^^s are Cauchy and so converge to F^JC. This guarantees that 
F^K is polyhomogeneous. The key observation here is that if V is a lift of V, then 

V (F, {K^^ - K'^)) = (V {K'' - K'')^ . 

The expression V {K''^ — K^.^) tends to because K'' ^ K as distributions conormal to 
H. □ 



Proof of Theoreml^ A simple wavefront set argument gives us that the K2f piece is smooth. 
Thinking functorially, we write 

In particular, if we write K = Ku, then 

wy'^ =u\ dxdy\ = {Pl)* ( -^(/3^^^)rg''^(ncf+ncf_)^''Vscfi^' 



We now use the decomposition K = Ki + K2 + , where Ki are as in Theorem HOI The 
lemma above shows that (5l is a b-fibration, and so we may use the pushforward theorem of 
Melrose in Lemma [9] to treat the contribution from K^. In particular, we note that if v has 
index set E at and vanishes to all orders at YL, then is polyhomogeneous on Xq, 
smooth at lcf+ and If, with index set E at ff+, and vanishing to all orders at ff_, lcf_, and 
scf because it is forward-directed. Index sets add when functions are multiplied, and so we 
know that rg''^(ricf+ricf_)^^^rscf-ft'3 ■ has index family Q given by 

Gff+ = F, Gq_ = 

Gif+ = {{s±{\) + m, 0) : m e No}, Gx,_ = 

Grf_ = 0, Grf+ = 

G'icf+ = {(j,0:^<J,jeNo}, Girf_=0 
Gscf = 0. 

Gif_^ must be modified when s+(A) — s_(A) is an integer. We may now use Lemma [9] to 
conclude that K^^f is polyhomogeneous on X with index set 



F U Gif+ U G 



Icfj 



To handle the contribution from ^^2, we use the lemmas above. Indeed, Lemmas |42| |44 
and US show that F'2/ is polyhomogeneous with index set Gicf+ U F. 

This leaves only the contribution from Ki. This is just a consequence of the local theory 
of paired Lagrangian distributions. Indeed, the work of Joshi in [llj implies that the push- 
forward of KiP^f exists as a smooth function. The uniformity of Ki down to the front face 
ff then tells us that Kif is polyhomogeneous with index set F. □ 
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18. An U' ESTIMATE 

As another apphcation of Theorem [381 we consider the behavior of the norms of a 
family of smooth compactly supported functions with support tending towards Y^. 

Suppose first that R e ^o°°(X) is a smoothing 0-pseudodifferential operator in the 
small calculus of [H] or (TH], supported near ff_|.. In other words, the Schwartz kernel of 
i? is a smooth function on Xq, supported away from If and rf and near ff_|_. Concretely, 
let is a smooth, compactly supported function on M", supported near (1,0), and with 
/o°° /r"-i ~ ^- Let X be a smooth function on X that is identically zero near F_. If 

R{s, z, X, y) = z)x{x, y), then R G \&o °°(^) is such an operator. 

Define now a family of compactly supported functions / given by 



f{x,y){x,y) = rI^, ^^-z-^,x,y] = Rv, 
\x X J 



where v = 6{x' — x)6{y' — y). Each f(x^y) is a smooth function on X supported in a compact 
neighborhood of {x,y) with unit L^{X; dg)-'n.oTm. 

Because f(x,y) is given by applying R to a. 6 function, pointwise bounds for E+f(^x,y) are 
equivalent to pointwise bounds on the Schwartz kernel of E^R. 

The following lemma is useful for obtaining pointwise bounds. 

Lemma 46. Suppose that K = Ku, where K is a polyhomogeneous function on Xg with 
index family T and supported near ff+. Suppose that J-" satisfies 

Fs+ = {(Sff + : J e No,p < Pj}, 
Fm+ = {(sicf + : j e No,p < Pj}, 
i^if+ = {(sif + J,p) : J e No,p < Pj}, 
Fri+ = {(srf + : j e No,p < Pj}, 

andpo = for each index set (i.e., no log terms appear in the top order part of the expansion) . 
Suppose further that 

n 1 

Sff > -, Sicf > -, Sif > 0, Srf > 0. 



Then K (considered as an operator) satisfies 

\\^'^\\l°^{x) — 

If > — f instead, then K is hounded L^{X; dg) L°°(X). 



Proof. The half-density Idgidg^il^^'^ lifts to a nonvanishing smooth multiple of 'r\I^^T^'^v 
near ff_|_, so we may write 



kv = kr;^l'\Z"^* (l d9LdgR\"^ 



where (3* is the puUback by the blow-down map Xq — > X x X. The assumptions above mean 
that kr^^^'^r^^^'^ is a bounded function on Xq and so is the puUback of a bounded function 
on X X xt 

The second statement of the lemma follows from the observation that dg = x^ dg. □ 
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Let Ki, K2, and K3 be the decomposition of given in Theorem [381 The Schwartz kernel 
of KjR is given by 

(27) [ K. ( — , ^ - + ^ + ^(5' / 5 ^) ds' dz' 

J \s'x s'x J 

where s' = ^ ^' = ^T^- This corresponds to writing the composition of operators A 
and 5 on X X X as 

KAB{x,y,x,y) = / KA{x,y,x',y')KB{x',y',x,y)dx' dy', 



where ka here denotes the Schwartz kernel of A. If ka and kb are instead functions of 
s, z, X, y, where s = | and z = this becomes 

/ X y-y' , A fx' y' - y ^ ~\ , / , / 

\x' x' J \X X J 

fx y-{y + xz') ,^ \ , , ^ , , 

Ka — r, ; , s x,y + XZ ] Kb is , z ,x,y) ds dz , 

\ s'x s'x I 



which yields equation (I2j 

Note that K\,K2, and i^s all vanish identically near rf+. Because the family j{x,y) is 
supported away from y_, we may assume that Ki are supported away from the "minus" 
faces ff _, lcf_, rf_, and scf. 

Consider first Kx, the piece corresponding to the paired Lagrangian singularity of E^. 
By Theorem [381 Schwartz kernel of K\ is given by {x'^~'^l'^K\, where K\ is a paired 
Lagrangian distribution on Xq supported away from If and rf . 

Because the fibers of integration in equation (127|) are transverse to the diagonal and to the 
light cone, it follows that K\R = where ui is a smooth function on Xq supported 

away from If and rf. In particular, ui is bounded, and Lemma [16] implies that 

||^l/(x,y)||ioo < C. 

Because R is defined only near ff+, we are free to use the improved parametrix of Theo- 
rem [38] rather than distribution of Theorem [lO] In particular, the symbol of the conormal 

1/2 

distribution K2 is bounded by r^J^ and so K2R satisfies the conditions of the second part of 



Lemma [IH] and so 

||^2/(i,y)||^oo(jf) < C- 

_g_ 

Consider finally the polyhomogeneous term K3. The Schwartz kernel of K^R sat- 

isfies the conditions of the lemma, and so 

||-^3/(x,y)||^s_(A)^^(^) < C. 

Now let / = max (0, — 3f?s_(A)). Putting the estimates for Ki, K2, and K3 together yields 

ll^+/{*.S)|L-iL°°{X) - 

An estimate for this family follows from Theorem [6] (from [23]). By reversing the roles 
of y_ and Yj^ in this theorem, we conclude that for any forward-directed / G x~^L^(X) and 
r > max (i, /(A)), there is a unique u G x~'^Hq(X) such that P{\)u = /, and 

\W\\x-^L2{X;dg) — \\'^\\x~'~H^{X) — ^ W f \\x-^ L2{X; dg) ■ 
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Here Hq is the 0-Sobolev space of order one, i.e., it measures regularity with respect to the 
xdx and xdy vector fields. 

In order to apply this estimate to the family f\x,y), it is important to understand how the 
L^-norms of the functions vary. Indeed, a simple calculation shows that 

2 f dx dy f ds dz 

where Cr depends on r, but not on x or y. In particular, WlW^-ri'^^x-dg) ~ ^''^^'^Cr- 
We may now prove Theorem [3l 

Proof. Interpolating between the L°° and estimates provides an estimate for p G (2, oo). 

Indeed, if r > max (^i, + , / = max (^0, + + , and i 

e e [0, 1], then 

(28) ||i?+/(£^y)||^_^a_,(i_9)^2/e(X;dg) — ' 

which finishes the proof. □ 

Remark 47. The proof of Theorem [3] uses the inclusion C L^. In particular, we ignore one 
derivative of and so we could modify equation (1281) to include a fractional derivative. 
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